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We discovered uniaxial oriented centers in silicon carbide having unusual performance. Here 
we demonstrate that the family of silicon-vacancy related centers with S= 3/2 in rhombic 15R- 
SiC crystalline matrix possess unique characteristics such as optical spin alignment existing at 
temperatures up to 250°C. Thus the concept of optically addressable silicon vacancy related centers 
with half integer ground spin state is extended to the wide class of SiC rhombic polytypes. The 
structure of these centers, which is a fundamental problem for quantum applications, has been 
established using high frequency ENDOR. It has been shown that a family of silicon-vacancy related 
centers is a negatively charged silicon vacancy in the paramagnetic state with the spin S= 3/2, VsG, 
perturbed by neutral carbon vacancy in non-paramagnetic state, Vc°, having no covalent bond with 
the silicon vacancy and located adjacently to the silicon vacancy on the c crystal axis. 


Spin centers in Silicon Carbide (SiC) have recently 
been put forward as favorable candidates for a new gen¬ 
eration quantum spintronics and sensorics, as well as 
quantum information processing because of the unique 
properties of their electron spins which can be optically 
polarized and read out [l-ll|. Particularly the optical 
control of the single defect spin in SiC has been real¬ 
ized at room temperature for the first time on the well 
studied silicon vacancy-related center with ground spin 
state S= 3/2 in 4H-SiC [loj]. The centers of the same 
origin persist also in 6H-SiC polytype and were proposed 
for use in quantum magnetometry [6( and room temper¬ 
ature operated MASERs [§|. Thus the family of the V Si 
related centers can be considered to form the basis of 
the SiC-based quantum sensors and quantum spintron¬ 
ics. The first solid state system on which optical detec¬ 
tion and manipulation of the single spin has been shown 
is the negatively charged nitrogen-vacancy (NV _ ) cen¬ 
ter in diamond. There are two main reasons why suc¬ 
cessful realization of single spin control in this system 
was possible. First, the NV _ centers possess high op¬ 
tically detected magnetic resonance (ODMR) contrast 
and their spin alignment persists at elevated tempera¬ 
tures 0, thusthe polarized spin state can be readout 
with high fidelity even on the single spin jl3| . Second, 
the precise atomic structure of the center has been elu¬ 
cidated Q0. This gave rise to the precise technol¬ 


ogy of NV’s production in predetermined topology and 
concentrations 0,0, allowing one to organize precise 
control of the NV’s electron spin interactions with the en¬ 
vironment. This permits to develop NV~ based quantum 
registers 0, [19] and quantum sensors [2C|], f2lj . 

Therefore for the successful development of the SiC 
two major issues have to be addressed. On the one 
hand spin centers with high breakdown characteristics 
and high ODMR contrast must be explored, on the other 
hand the proper model for these spin centers must be es¬ 


tablished. We focused our efforts on these two tasks and 
present our results in this work. 

SiC single crystals of 15R polytype, were grown by the 
Physical Vapor Transport method. The concentration of 
uncompensated nitrogen was in the range «10 16 cm -3 . 
Vacancies were introduced by irradiation of the crystal 
with 1.4 MeV electrons with a dose of 10 18 cm -2 . The 
crystal structure was confirmed by means of Raman scat¬ 
tering spectroscopy 0 , HH ■ Results are summarized in 
Supplementary Information (SI). 
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FIG. 1. (Color online) EPR spectra of the 15R-SiC taken 
at temperatures T= 303K, T= 479K, T= 523K. Insets show 
light-induced inverse population of the spin sublevels of V2,V3 
and V4. 

First we measured a cw X-band (a: 9.3GHz) elec¬ 
tron paramagnetic resonance (EPR). Fig.l shows spectra 
recorded at three different temperatures (up to 250° C) 
under continuous laser light illumination with A= 785nm 
and the magnetic field oriented parallel to the c - axis of 
the crystal (B11 c) (without optical excitation the signals 
at room temperature are near the noise threshold). The 
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TABLE I. The values of D, g , wavelength of ZPLs and Ti 
and T 2 ESE measured at room temperature 


center 

D(MHz) 

g 

ZPL (nm) 

Ti fis 

T 2 g. s 

V2 

69.6±1 

2.005(1) 

886.5 

80 

10 

V3 

5.8±1 

2.005(3) 

904 

120 

10 

V4 

25.1±1 

2.005(3) 

917 

105 

9 


spectra consist of three pairs of EPR transitions V2, V3, 
V4 attributed so from the corresponding Zero Phonon 
luminescence Lines (ZPL) as will be shown further. The 
angular variation of the splitting value between each pair 
of lines (see SI) can be fitted by the axially symmetric 
spin Hamiltonian with parameters summarized in Table 
I: 

H = 9 hbB 0 -S + D(S 2 z -^S(S+1)) (1) 

where the first and the second terms correspond to the 
Zeeman interaction and fine structure splitting, respec¬ 
tively, Hb and g is the Bohr magneton and electron g 
factor, respectively, total electron spin is S= 3/2 and D 
is the fine structure parameter, e.g., zero-field splitting 
(ZFS) = 2D. The evidence of the S= 3/2 state will be 
provided from our ODMR measurements. Thus, all three 
centers are c-axis oriented with C 3 „ point-group symme¬ 
try. 

A phase reversal was observed for one of the two tran¬ 
sitions in each pair of lines. Such behavior in the emis¬ 
sion/absorption mode of the microwave power can be ex¬ 
plained by optically induced spin alignment of the spin 
sublevels of the centers. Note the unique properties of 
these centers, where optical polarization can be created 
at high temperatures up to (T« 250°C), which indicates 
that the spin system of the centers is well isolated from 
the crystalline lattice vibration. In this aspect they are 
similar to the NV - centers in diamond, with important 
exception - the ZFS of each center stays stable at the 
whole temperature range used while the NV-centers ex¬ 
perience the huge thermal shift H3,0,0. 

To establish the model of the optical polarization of 
the V2, V3 and V4 centers we determined the ordering 
of their spin sublevels (the sign of the ZFS). Low temper¬ 
ature EPR measurements, conducted at high frequency 
(95GHz), and Electron Nuclear Double Resonance (EN- 
DOR) study allowed us to determine that ZFS is positive 
for each center (See SI) and establish the scheme of the 
optically polarized spin sublevels as shown in the insets 
of Fig.l. 

We then decided to probe the coherence times of the 
ensembles of the spin centers using EPR X-band pulsed 
technique. The results are summarized in Table I for 
each center. The spin lattice relaxation times (Ti) and 
spin-spin relaxation times (T 2 ) (see Table 1) turned out 
to be of the same order of magnitude with those of 4H- 
SiC spin centers measured previously on the single silicon 


vacancy related centers 0 - 




FIG. 2. (Color online) (a) PL spectrum of the 15R-SiC. 
ZPLs are labeled as VI, V2, V3, V4. (b) cw X-band DD- 
EPR spectra obtained at temperature T= 10K and (B||c). 

We then proceeded to measure optical properties of the 
centers discussed above in order to obtain the resonant 
optical control of their spin states. For such investiga¬ 
tions the low temperature photoluminescence (PL) spec¬ 
trum was measured (Fig.2 (a)). To find the zero photon 
line (ZPL) that corresponds to each of the EPR centers 
shown on the Fig. 1 we excited each ZPL with laser flash. 
Excitation of each ZPL gives rise to one pair of EPR lines 
with ZFS corresponding to those determined from our 
cw-EPR measurements (Table 1). Resulting spectra are 
summarised in Fig.2 (b) and Table 1. Particularly ex¬ 
citation into V4 line with energy E= 1.3522 eV induced 
the signal with ZFS equal to 50 MHz, in V3 (E= 1.3716 
eV) and V2 (E= 1.3989 eV) with ZFS 13 MHz and 138 
MHz, respectively. So we have shown that each center in 
15R-SiC has its own optical fingerprint. Excitation into 
VI line didnt give rise to any EPR signal. 

For any application in spin center-based spintronics 
or quantum sensorics one have to get access to effective 
readout of polarized spin state. We further demonstrate 
that V2, V3, V4 centers can be optically addressed and 
readout with high fidelity at room temperature by means 
of a standard ODMR technique. 

The ODMR spectrum in the Figure 3 shows relative 
change of the photoluminescence intensity APL/PL as 
a function of applied RF. Resonances at v = 13.1 MHz, 


























3 



Frequency (MHz) 


FIG. 3. (Color online) Lower spectrum shows ODMR signal 
of the V2, V3, V4 detected under laser excitation A= 785 nm 
in 0.05 mT external magnetic field applied in order to com¬ 
pensate the influence of Earth magnetic field. Experimental 
frequency dependence of the ODMR signals as a function of 
magnetic field are shown for each center. Dashed lines are 
the calculated dependence for the allowed transitions using 
equation (1), S= 3/2 


v = 50.6 MHz and u = 138 MHz were observed that agree 
well with ZFS of the centers V3, V4 and V2 determined 
from our EPR. measurements. 

Using spin Hamiltonian from the Eq. (1) we further 
calculated that the ground states of the V2, V3 and V4 
have a spin state S= 3/2. Calculated frequency depen¬ 
dencies of the ODMR signals at different magnetic fields 
(B11 c) coincided well with observed spectra (Figure 3). 

The unusual properties of spin centers in 15R-SiC ob¬ 
served in this work are several times higher compared 
to similar signals in the 6H- and 4-H polytypes. Tak¬ 
ing into account high temperatures at which the optical 
polarization persisted (up to 250°C) and strong ODMR 
contrast these centers seem to be favorable candidates 
for the spin center-based quantum applications. More¬ 
over, changing the excitation energy and the optical reg¬ 
istration window, the ODMR contrast can be further in¬ 
creased. When ODMR signal was measured in 15R-SiC 
under excitation of the V4 center with laser A= 808 nm 
we observed about five times more intense signals than 
in other polytypes (unpublished results). 

From the experiments described above it becomes clear 
that the V Si related centers with optically addressable 
S= 3/2 ground states exist not only in hexagonal but 
also in rhombic polytypes. In order to advance under¬ 
standing functionality of the observed spin systems we 
determined the microscopic model of the centers. Cur¬ 
rently there are two such models under consideration. 


The first model postulated these centers to be the low 
symmetry configuration of the V Si _ with S= 3/2 (V Si _ ) 
[HI, 0]. According to the second model these centers 
are the V Si “ S= 3/2, perturbed by neutral carbon vacan¬ 
cies in nonparamagnetic state, V c °, located adjacently to 
the silicon vacancy V Si “ on the c - axis of the crystal, see 
Q and references therein. 

In this work we present the strong evidence in favor of 
the second model by means of Electron Spin Echo (ESE) 
detected high frequency ENDOR. 

To this end we measured the W-band (95 GHz) EPR 
spectra of the V2, V3 and V4 centers under A= 785 nm 
laser excitation for two orientations ( 9) of the magnetic 
field B with the respect to the c axis (Fig.4 (a)). We 
here focused on the V2 center with larger ZFS= 138 MHz 
to obtain its ENDOR spectra. Allowed dipole magnetic 
transitions with A M$= ±1: 3/2 o 1/2 are indicated by 
abbreviation If (low field), —3/2 o —1/2 indicated by 
hf (high field). The ENDOR spectra were obtained by 
monitoring the intensity of the ESE, following three mi¬ 
crowave 7 t/ 2 pulses, as a frequency function of the pulse, 
applied between the second and third microwave pulses 
[30( (Fig.4 (b,c)). To describe the liyperfine (HF) inter¬ 
action with nuclei of silicon and carbon, which reside in 
different shells, one should add to the Spin Hamiltonian 
given in Eq. (1) the term: 

ESAJ, + ESAjlj + ESA/.I k, where A, and A/. are 
tensors which describe the HF interaction with i-th Si and 
k-th C atoms located in different neighbor shells of the 
Si sites, Aj is tensor that describes the HF interaction 
with j-th Si atoms located in neighbor shell to the C site. 

The ENDOR transition frequencies determined by the 
selection rules AMs=0 and AM/ =±1 are given by [3lJ 

vendori = h~ 1 \Ms(a i + bi(cos 2 9 — 1)) - g n i^ n B\ (2) 

Here at and bi denote the isotropic and anisotropic part 
of the interaction with the «-th nucleus (e.g., for Si), re¬ 
spectively, 9 is the angle between the direction of the 
external magnetic field B and the hyperfine interaction 
tensor, g n i^niB/h is the nuclear Zeeman term or the Lar- 
mor frequency f l, g n i and g. n i are g factor of nucleus i 
and its Bohr magneton, respectively ( g n is negative for 
29 Si and positive for 13 C). For axial symmetry HF inter¬ 
action in terms of the principal values are given by A|| = 
a + b, Aj_ = a - b. 

Our data suggest the structural model of the V2 cen¬ 
ters as shown in Fig.5. The model represents two possible 
configurations (V Si - - V c °) and (V Si _ - V c °) . First, 
prior to the analysis of obtained ENDOR data, we would 
like to discuss the HF interactions directly observable 
from the ESE spectra in Figure 4 (a). 

The hyperfine structure,shown in the enlarged scale on 
the left inset of Fig.5(a), arises from the HF interactions 
with 29 Si in the next nearest neighbor (NNN) shell to 
the V Si - site, and is characterized by the constant A= 
2.97G. The similar interaction was observed by others in 
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FIG. 4. (a) (Color online) Optical induced ESE spectra ob¬ 

tained in the 15R-SiC at two angles ( 9 ): 9 = 0° and 9 = 12,5°. 
(b), (c) ENDOR spectra measured at the If and lif transitions 
indicated in (a). 29 Sii and 29 Si 2 - 4 indicate ENDOR lines cor¬ 
respond to the presence of 29 Si atoms in the NN shell with 
respect to Vc°- One and two asterisks in (b) show transi¬ 
tions corresponding to ft + 1 /21Ay | and f l + 3/2|Aj| for HF 
interactions with axial NN Si atom with respect to Vc°- 


4H- and 6H-SiC [27], |32j. Note that changing the ori¬ 
entation of the magnetic field does not induce the line 
shifts; meanwhile the strong anisotropy of the ESE line 
width was observed. On the other hand, the HF inter¬ 
actions with nucleus located in the first NN shell to the 
V Si “ site are known to be strongly anisotropic and re¬ 
flect the tetrahedral symmetry of the nuclear spin loca¬ 
tions. Hyperfine lines arising from interactions with 13 C 
are shown in Fig.4 (a). 13 Ci denotes interaction with the 
carbon atom oriented along the c - axis and the 13 C 2-4 
denotes the interactions with carbon atoms located in 
the basal plane with bonds inclined by the angle 9= 71° 
relatively to the c - axis. These carbon atoms are shown 
schematically in Fig.5. The HF structure arising from 
such interactions is characterized by the parameters A| 
= 30.2 G (84.6 MHz) and Aj_ = 12 G (33.6 MHz) which 


match closely previously reported values for the V Si _ cen¬ 
ters in 4H- and 6H-SiC [ 23 , fril - frit Therefore our ESE 
measurements point to strongly anisotropic nature of HF 
interactions of the V Si _ electron spin with carbon nuclear 
spins in the NN shell and reflect their tetrahedral config¬ 
uration. 

From the measured ENDOR spectra (Fig.4 (b,c)) all 
the lines observed can be attributed to the HF interac¬ 
tions of V Si “ with 29 Si. To reveal this we need to con¬ 
sider Eq.2, which predicts that for S= 3/2 each nuclei 
i induces two sets of ENDOR transitions located at a 
distance 3/2A* and 1/2A, for the If EPR line, -3/2A i 
and -1/2A,; for the hf, from the corresponding Larmor 
frequency f l- Thus, the ENDOR lines observed are due 
to HF interactions with 29 Si, because the fj, of the 29 Si 
(1= 1/2, abundance 4.7%) is about 28.2MHz. The more 
accurate value of the U depends on the magnitude of 
the magnetic field for the EPR transition and is marked 
by arrows in Fig. 4(b,c). HF interactions of the V Si _ 
electron spin with 29 Si nuclear magnetic moments in the 
NNN shell observed on the ESE spectra can also be seen 
in our ENDOR spectra (Fig. 4 (c), labeled 29 Si(NNN)). 
Position of these lines corresponds to the positive den¬ 
sity of the electronic wave function around the silicon va¬ 
cancy. Unexpectedly we also observed the lines that were 
characterized by the negative spin density on Si nucleus. 
In Fig. 4(b) we show the expanded part of the ENDOR 
spectra for lines that correspond to the negative spin den¬ 
sity measured at two different orientations of B relative 
to the c - axis (9= 0° and 9= 12.5°). Strong anisotropy 
of the signals (shown by shadowed area in Fig.4(b)) was 
observed. Earlier in discussion of ESE spectra we showed 
that interactions with 29 Si located in NNN and more dis¬ 
tant shells are almost isotropic, meanwhile, interactions 
with 13 C located in the NN shell possess anisotropy. In 
ENDOR spectra we observed for HF interactions with 
with negative spin density 29 Si the same anisotropic de¬ 
pendence as that for 13 C located in the NN shell seen in 
ESE. To explain such an anisotropy, one has to identify in 
the model (Fig. 5) Si atoms which reflect the same sym¬ 
metry as that for C atoms in NN shell of silicon vacancy. 
For Si atoms, such a configuration could be found only 
at the tetrahedron vertices around the carbon site in the 
SiC lattice. Position of ENDOR lines labeled as 29 Sii and 
29 Si 2-4 in Figure 4(b) agrees well with proposed config¬ 
uration and reflects the HF interactions with axial (Sii) 
and basal (Si 2 - 4 ) nuclear spins. The line correspond¬ 
ing to HF interactions with axially oriented Si nuclei for 
transitions with M$= 3/2 is shown on Fig.4 (c) marked 
by two asterisks. The constants of HF interactions with 
negative spin density are relatively large and correspond 
approximately to the A|| = 2.2 MHz (ss0.8 G) Aj_ =1.3 
MHz (/« 0.5 G). Notably these constants describe well 
the anisotropy of the linewidth observed on ESE. 

HF interactions with negative spin density 29 Si could 
be explained if the spin density was located on four Si 






































5 


nuclei placed around the carbon vacancy being in the 
non-paramagnetic neutral charge state, V c °. Therefore 
we conclude that both the V c ° and paramagnetic S= 
3/2 V Si “ form the V2 center. Negative spin density is 
caused by the spin polarization (similar to the core po¬ 
larization for transition metals (3a . ®). This spin polar¬ 
isation arises from exchange interaction with S= 3/2 that 
leads to partial decoupling of coupled covalent bonds in 
the V c ° site. The other interactions with smaller nega¬ 
tive spin density observed on ENDOR spectra could be 
attributed to HF interaction with more distant Si nu¬ 
clei from the V c ° site. V c ° distorts the crystal lattice 
which in turn lowers symmetry of the V Si _ , so the car¬ 
bon vacancy could be expected at the center with the 
largest ZFS (V2) in the position most close to the V^T. 
It is important to note that the double number of lines 
observed in the ENDOR spectra points to the presence 
of two similar centers with slightly different parameters 
of HF interactions (Fig.4 b, 8 = 12.5°). In 15R-SiC it 
could be connected with the unit cell structure. Thus, the 
model of the V2 center could be reconstructed as shown 
in the Fig.5. The ENDOR spectra of (V3) and (V4) in 
15R-SiC (not shown here) exhibited the same pattern of 
the HF interactions with 29 Si nuclei, having slightly dif¬ 
ferent values of hyperfine splitting, so the structure of 
these centers has to be similar to that of V2 (V Si “ 
V c °). The only difference would be that the distance be¬ 
tween carbon-silicon vacancies along the c - axis should 
be larger than that of V2. Moreover, we have shown 
here that the properties of V Si _ related centers with S= 
3/2 in 15R SiC resemble those of 4H, 6H , therefore the 
model proposed here for V2 in 15R (Fig. 5) could be 
extended for hexagonal polytypes. Supporting evidence 
for possible model universality will be provided in later 
publications. 



FIG. 5. (Color online) The fragment of the 15R-SiC struc¬ 
ture with two hexgonal (hi, h2) and three quasicubic (kl, k2 
and k3) sites for silicon. Vacancies are labeled by squares. 
Carbon and Silicon atoms in the NN shell to the Vsi - and 
V c ° are labeled as Ci _4 and Sii_ 4 , respectively. Structure 
of the V2 center is indicated as (Vs A - Vc°) and (Vsi - 
Vc°)". 

In conclusion in this work we have provided evidence 
for the presence of three uniaxial spin centers in 15R-SiC 
lattice. These centers with distinct ZPL and ZFS were 


characterized by long spin coherence times at room tem¬ 
perature. Their optical spin alignment persisted up to 
250°C and the value of the ZFS did not show temperature 
dependence. ENDOR study revealed the presence of two 
types of HF interactions of the electron spin of these cen¬ 
ters with surrounding 29 Si nuclear spins, characterized by 
negative and positive spin densities. This allowed us to 
establish the model for the V2 center. This knowledge 
is important for understanding the processes leading to 
the decoherence [37] and paves the way to the effective 
and precise production of silicon vacancy related centers. 
The newly discovered silicon vacancy-related centers are 
well isolated from the SiC crystal matrix, therefore these 
centers lend themselves for further development of SiC 
based quantum applications and quantum registers based 
on coupled electron spins. Furthermore, our results ex¬ 
tend the concept of optically addressable V Si _ related 
centers to the rhombic polytypes of silicon carbide. 

This work has been supported by MES of Russia 
No 14.604.21.0083 and by the RFBR No 13-02-00821. 
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